Azimuthal decorrelations of dijets in QCD 
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We report on the status of the QCD analysis of dijet azimuthal decorrelations. We 
emphasise the relevance of resummation of soft and collinear enhancements in describ- 
ing these observables in the region where the two jets are nearly back-to-back in the 
transverse plane. We also discuss the sources of theoretical uncertainties and possible 
research directions aimed at their reduction. 

1 Dijet azimuthal decorrelations 

The study of jet cross sections constitutes one of the most important investigation tools 
for the physics of hadronic final states. These observables are favoured in many aspects, 
both theoretically and experimentally. First they are infrared and collinear safe|f] which 
makes it possible to compute them at all orders in QCD. Hadronisation corrections are 
quite small, being suppressed by inverse powers of the transverse energy of the jets, and their 
dependence on the jet radius R is under control [3J. Methods to eliminate the contamination 
of underlying event and pile- up are also being developed [3J. Furthermore, measurements 
involving jets show less experimental uncertainties with respect to the corresponding ones 
involving particles (see for instance Ref. [I] for a recent analysis). 

Among jet observables we consider the azimuthal decorrelation of a pair of jets [5l [3]. 
At tree level two jets produced in a hard collision are highly correlated, their transverse 
momenta with respect to the beam being exactly back-to-back, which implies that their 
relative azimuthal angle is equal to tt. Additional QCD radiation decorrelates the dijet 
system, moving away from its Born value. This feature, together with the fact that any 
measurement involving angles, i.e. ratios of momenta, is less sensitive to jet energy scale 
uncertainties, make azimuthal decorrelations an ideal ground to explore all-order properties 
of QCD dynamics. In particular, in the region where A</> is close to tt, we expect the 
dominance of configurations with multiple soft and collinear emissions. On the contrary, 
when A(f> is small, if the available rapidity range is large enough, one may expect the 
occurrence of multiple hard emissions decorrelating the two jets. Besides these perturbativc 
features, azimuthal decorrelations are also affected by intrinsic transverse momentum of 
partons inside the incoming hadrons, thus making it possible to explore this non-perturbative 
aspect of QCD. 

As far as the experimental situation is concerned, dijet azimuthal decorrelations have 
been measured both at the Tevatron [5] and at HERA [6j. In hadronic collisions, one observes 
that next-to-leading order (NLO) QCD predictions obtained with NLOJET++ [7] fail for 
Acj) near tt. This region is instead described quite well by the Monte Carlo event generators 
HERWIG [H] and PYTHIA [5]. This reveals that at the Tevatron, where the momentum 
fraction x of struck partons is quite large, multiple soft and collinear emissions are able to 
account for the azimuthal decorrelations in the quasi-back-to-back-region. Furthermore, due 



a Some infrared unsafe cone algorithms are still on the market, but efforts are being made to move to 
infrared safe jet definitions [TJ. 
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to the sensitivity of Monte Carlo predictions to changes in the shower evolution parameters, 
these observables should be exploited for the tuning of event generators [5] . HERA instead [5] 
explored the region of small x, where large deviations from Monte Carlo predictions have 
been observed. However, one finds better agreement with NLO, except of course in the 
region A</> ~ tt. A question then arises on whether the Monte Carlo tuning performed at the 
Tevatron at large x will be able to describe LHC data which involve smaller values of x. An 
attempt to predict azimuthal dccorrclations both at large and small A(j> is represented by 
CASCADE [10], an event generator producing hard QCD emissions according to the CCFM 
equation [TTJ. CASCADE, as expected, describes the azimuthal decorrelation in the region 
of small A0, while fixing the distribution at A(j) close to tt by means of an unintegrated 
parton distribution containing information on the transverse momentum of the incoming 
parton inside the proton. 

We decided to take another point of view. First we observe that a resummation of soft 
and collinear emissions gives rise to a distribution that rises to a constant value for A(f>— >7r, 
consistent with what is seen in the data [3 [6]. Moreover, the agreement of NLO QCD with 
data in the region of small A<p suggests that in the x region probed at HERA only few extra 
hard emissions could be considered. It is then reasonable to investigate whether HERA data 
can be described by a resummation of soft/collinear enhancements matched to exact NLO. 

2 Resummation of soft and collinear logarithms 

The resummation of soft and collinear logarithms is performed by considering the observable 
S(A), the probability that \tt — A(j>\ < A. The azimuthal decorrelation can be then obtained 
by differentiating S(A). The rate E(A) contains logarithmic contributions up to a™ In 2 ™ A, 
which become large in the back-to-back region A <C 1, and have to be resummed at all orders. 
The first step to perform such a resummation is to investigate the behaviour of [tt — A(f>\ 
after a single soft emission fc, collinear to each of the hard legs (one or two incoming and 
two outgoing). If k is collinear to any of the incoming partons, we have 

kt 

I vr — A(j)\ ~ 1 sin^l , 

Pt,i 

where k t is the emission's transverse momentum (with respect to the beam in hadronic 
collisions or the virtual photon momentum in DIS), </> its azimuthal angle, and pt,i is the 
transverse momentum of the highest-pt jet. This result does not depend on either the jet 
algorithm or the recombination scheme of particles within a jet. If instead one considers an 
emission collinear to one of the outgoing legs, in any jet algorithm, soft gluon k and its parent 
parton will be clustered in the same jet, and the result will depend on the recombination 
scheme. In any scheme that adds three-momenta vectorially (i?-scheme, P-schemc, E0- 
scheme) the azimuth of the jet does not change, so that \tt — A(j>\ is completely unaffected by 
emissions inside the hard jets. Therefore the azimuthal decorrelation turns out to be a non- 
global observable [12] , its resummation becomes extremely tricky and will not be discussed 
here. If instead one uses a pt-weighted recombination scheme, as is done at HERA, we have 

kt kt 
Itt- A(f>\ ~ — Isin0-0|, iejetl; \ir - Ac/)\ ~ — | sin0 - (tt - 0)1 , i £ jet2. 
Pt,i Pt,i 
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For an arbitrary set of secondary soft and collinear emissions, we have |13| 



A0| 
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z^jets 
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i£jetl 



Pt,l 
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and the observable is now global. We stress that the choice of the pt-weighted scheme is 
observable specific. For other observables a different choice could be needed in order to 
achieve the best theoretical accuracy. 

Resummation is better performed in impact parameter space, where the result for E(A) 
schematically reads [T3] 



2 r°°r]k 
E(A) = - / — sin(6A) 



/a(/W b ) 



n 



t(b) 



5(6). 



A = A, 



(1) 



Each function in the above equation accommodates all-order real or virtual corrections in 
selected phase space regions. The ratios of parton densities f a {jJLF /£>)/ /o(^f) embody real 
and virtual contributions up to the scale ^_f/&, above which real emissions are forbidden 
by the observable definition. All remaining virtual corrections are included in the radiator 
functions R m {b) 1 R ou t(b), containing virtual corrections collinear to incoming and outgoing 
legs respectively, and the soft function S(b), accounting for soft large-angle gluons. The 
above resummation is valid within next-to-leading logarithmic (NLL) accuracy, i.e. accounts 
for all terms a™ ln" +1 & and a™ ln n 6 in the logarithm of the Fourier transform of S(A). At 
first order in a s the expressions for R- m (b) and R on t(b) read 



C in - In 2 6, 

7T 



■Rout (&) 



Cout O-s , 2 i 

In o . 



where C- lri and C ou t are the total colour charges of incoming and outgoing legs respectively^] 
Note that R ln (b) is the usual radiator for px resummation (e.g. Drell-Yan lepton-pair pr [33]) 
while the expression for i? ou t (b) has never been encountered before and reflects the unusual 
dependence of \tt — A(f>\ on the rapidity and azimuth of each gluon emitted from outgoing 
legs. The soft function S(b) is in general a matrix in the colour space spanned by the hard 
emitting partons, and depends on emitters' four- velocities. 



3 Towards phenomenology 

We consider azimuthal dijet decorrelation in DI^| with the same kinematical cuts adopted by 
HI [B], that is two jets with —1 < r/ ab < 2.5 and transverse momenta in the hadronic ccntre- 
of-mass (HCM) frame p t > 5GeV. We give predictions for the differential cross section in the 
angular distance A(f> in the HCM frame. Results are shown in Fig. [1] for a selected value of 
Xb and Q 2 , reported in the figure. One can see that around A^> ~ 170° both LO and NLO 
predictions obtained from NLOJET++ diverge due to the presence of large logarithms. The 
NLL resummed curve in Eq. (fTJ) instead rises to a constant value, in qualitative agreement 
with what is seen in the data [BJ. 

b For instance in DIS with an incoming quark d n = Cf and C ou t = Cf+Ca- 

c This observable, denned with a pt-weighted recombination scheme, is global. This is not so for the 
corresponding observable considered at the Tevatron, which is non-global [5]. 
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Figure 1: Fixed order and resummcd predic- 
tions for azimuthal dijet dccorrclation in DIS. 



The resummed expression in Eq. JTJ) 
contains some ambiguities that have to be 
addressed before comparing to data. First 
all functions of b in the Fourier transform 
of E(A) result from integrals of the run- 
ning coupling a s (kt) down to kt = pt,i/b. 
When b becomes large, this value can ex- 
ceed the Landau pole, thus making Eq. JT]) 
ill-defined. We then decided to cut the b- 
integral at the value 6 max corresponding ex- 
actly to the Landau pole. However, since 
the 6-integral is strongly suppressed in the 
large-6 region, it is not affected significantly 
by reducing 6 max by a factor, as can be 
seen in Fig. [5] comparing the curves labelled 
b < frmax and b < 6 max /2. This means that 
even in the region A(p ~ 180° dynamics is 
not dominated by low momentum scales. At 
small b one has to take into account that set- 
ting 6 = corresponds to computing the total cross section. This imposes the requirement 
that the Fourier transform of E(A) tend to unity as b — ► 0. This is achieved by either 
freezing its value at one for b < 1 (this is what is done in the plot in Fig. [1]) or making 
the replacement b — ► yl + b 2 [H] ■ Figure [5] shows that either choice (corresponding to the 
curves b < 6 max and 6— >(1 + b 2 ) 1 / 2 ) produces almost indistinguishable results. 

The most important source of theoreti- 
cal errors is the treatment of f a (fJ-F/b), since 
the actual factorisation scale [ip/b can be 
smaller than the inverse size of the proton 
Rp 1 . The curve in Fig. [T] has been ob- 
tained by freezing the parton densities for 
b > fipRp. However, restricting the b in- 
tegral to b < [IfR p (see Fig. [2J gives a re- 
sult that can be 50% smaller in the large 
A(f> region. This suggests that to provide 
accurate predictions for azimuthal decorre- 
lations in DIS one may need to consider the 
effect of the intrinsic transverse momentum 
of partons inside the proton. 

Finally, in order to obtain quantitative 
predictions, one has to match NLL resum- 
mation to exact fixed order calculations. 
After matching, S(A) gets multiplied by 
a coefficient function (1 + C\a s + . . . ) in 
such a way that in its expansion all terms 

a™ ln 2n_2 A are correctly taken into account. The determination of C\ poses two theoretical 
problems. The first is that, in order to achieve such an accuracy, one has to be able to 
compute C\ for each hard colour configuration. This can be done in the soft/collinear limit 
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Figure 2: The impact of different prescriptions 
on azimuthal dijet dccorrclations in DIS. 
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via a flavour jet-algorithm [15], but requires that one modify NLOJET++ to include the 
information on the flavour of produced partons, as has been done in hadron-hadron collisions 
in Ref. [16]. This work is currently in progress. A second issue is the fact that, after sub- 
tracting all logarithms from the LO prediction, the result tends to a constant value which is 
large and negative, thus potentially giving a negative coefficient function. This reveals that 
probably part of the coefficient function will have to be resummed and eventually exponen- 
tiated, as has been done for the Drell-Yan total cross section [17]. Potentially large terms 
may contain 7r 2 , coming cither from Coulomb phases or from azimuthal integrations, or large 
logarithms of Bjorken x. The need for better understanding of the coefficient function is one 
of the reasons why we have decided to study these sub-leading contributions in the simpler 
case of the Z boson a<r distribution in hadron-hadron collisions [18j . The resummation for 
this observable is very similar to that for azimuthal decorrelations, with the simplification 
that it does not involve any outgoing jets and that it should have no contributions that 
are logarithmically enhanced by kinematics (e.g. by lnl/x). These last contributions could 
be addressed by extending soft collinear resummation to include also multiple hard gluon 
emissions as is done for instance in CASCADE. 
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